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ABSTRACT 

Two well-studied white dwarfs with helium-dominated atmospheres (DBs) each pos- 
sess less hydrogen than carried by a single average-mass comet. Plausibly, the wind rates 
from these stars are low enough that most accreted hydrogen remains with the star. 
If so, and presuming their nominal effective temperatures, then these DBs have been 
minimally impacted by interstellar comets during their 50 Myr cooling age; interstellar 
iceballs with radii between 10 m and 2 km contain less than 1% of all interstellar oxygen. 
This analysis suggests that most stars do not produce comets at the rate predicted by 
"optimistic" scenarios for the formation of the Oort cloud. 

Subject headings: interstellar medium - planetary systems - white dwarfs 

1. INTRODUCTION 

There are two reasons to imagine that interstellar comets could be widespread. First, they 
could be an important as-yet undetected reservoir of interstellar oxygen. Second, models for the 
formation of the solar system's Oort cloud predict that more comets are ejected into the interstellar 
medium than remain gravitationally bound to the Sun. If our solar system is at all typical, then 
the birth of stars and planets might commonly lead to the creation of many interstellar comets. 

Despite decades of effort, analysis of the rate of appearance of new comets in the solar sys- 
tem has not provided a strong constraint on the spac e density of interstella r comets. Substantial 
progress may require another approach. In this spirit, Shull &: Stern (j 19951 ) argued that one pos- 



sible explanation for the "low" observed frequency of gamma ray bursts arising from the Milky 
Way is that there are few interstellar comets impacting upon neutron stars. However, there are 
enough uncertainties in this approach, that their result that fewer than 4% of stars have Oort clouds 
similar to the Sun's is not completely firm. While also employing degenerate stars as a tool, here, 
we proceed somewhat differently. We argue that the small amounts of hydrogen accumulated by 
warm DB white dwarfs significantly limits the rate of iceball impact with these stars and therefore 
the space density of interstellar comets. Our method requires combining an understanding of the 
interstellar medium with the physics of external pollution of white dwarfs. We therefore subdivide 
this Introduction into subsections to describe the different threads of our approach. 
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1.1. Interstellar Oxygen 



Including both atoms in the gas phase and those estimated to be contained within dust 
grains, the abundance of oxygen within diffuse clouds, n(0)/n(H), is within 15% of 5.1 x 10~ 4 
(jCartledge et al. 1 12004 ) . While this value is in agreement wi th recent determinatio ns of the solar 



photospheric oxygen abundance o f n(0)/n(H) = 4.6 x 10 4 (jAsplund et al 



5.8 xlO" 4 (Caffau et al. 



20041 ) or n(0)/ra(H) 
2008), the tr ue solar abundance may be as high as n(0)/n(H) = 7.2 



x 10~ 4 ( Delahave Sz Pinsonneault 1 12006). Even with the problematical assumption that the in- 
terstellar and solar oxygen abundances are the same, the interstellar oxygen abundance is ~50% 
uncertai n. The possibility of a missing reservoir of interstellar oxygen has been considered for 
decades (jGreenberg Ill974l ). and here we re-assess the hypothesis that appreciable oxygen is car- 
ried within interstellar comets. Furthermore, although not directly related to the question of the 
total abundance of interstellar oxygen , the observed decrease of n(0)/ra(H) with column density 
(jCartledge et al. 1120041 ; I Jenkins 1120091 ) is not well understood an d presents a "c risis" in our under- 
standing of the behavior and distribution of interstellar oxygen (jWhittet Il2010i ). 



With an interstellar density of hydrogen nuclei near the Sun of 1.15 cm -3 ( Bohlin et al. Ill978l ) 
and with n(0)/n(H) = 5.1 x 10 -4 , the mass density of interstellar oxygen is 1.6 x 10~ 26 g cm -3 . 
Next, we compare this estimate with observational bounds and theoretical expectations for the 
space density interstellar comets and the amount of oxygen that they contain. 



1.2. Interstellar Comets? 



To date, all observed solar system co mets appear to have a lways been gravitationally bound 
to the Sun. Consequently IWhipple I (119751 ) and ISekanina I (jl976l ) placed upper limits of 2 x 10~ 26 



g cm 3 and 4 x 10 26 g cm 3 , respectively, t o the mass den sity of interstellar comets. Using the 



same approach but with updated parameters, iFrancis I (120051 ) placed an upper bound to the space 
density of int erstellar comet s of 4.5 x 10 -4 AU -3 . With an average mass per comet as much as 
1.2 x 10 18 g (IFrancis 1120051 ). and assuming about 50% of a comet's mass is oxygen (jGreenberg 



1998), then the upper bound to the mass density of oxygen within interstellar comets is 8 x 10 



-26 



g cm 3 - larger than the amount of gas-phase interstellar oxygen ( Stern fc Shull Ill990i ). 



Conventi onal dynamic models for the formation of the Oort cloud from objects in the inner 
solar system (jDuncan et al. 1119871 ) predict that there could be as many as 50 comets ejected into 
the interstellar medium for every com et that remains bound to the Sun. Unconventional models for 



the Oort cloud (jLevison et al. Il2010i ) also must require significant numbers of interstellar comets. 



An "optimistic" model for the space density of interstellar comets presumes that all stars form 
massive comet clouds similar to the solar system's Port cloud an d that most of these comets are 
lost into the interstellar medium. With this picture, IStern I (jl990i ) predicted that the space density 



of interstellar comets could be as high a s 2 x 10 cm 3 (or 7 x 10 AU ), comparable to 



the observational upper limit of IFrancis I (j2005l ). Consequently, "optimistic" models predict that 
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occasional interstellar co mets can be observed within the solar system ([McGlvnn fc Chapman 



1989 : ISen &: Ramal ll993). However, since the mass of the Oort cloud, the fraction of comets that 



are ejected into the interstellar medium and t he fraction of stars that create analogs to the Port 
cloud are each unc ertain by a factor of ~10 (jStern Ill990l ; IShull Sz Stern 1 1 19951 : IWeissman 1 1 19961 : 



Pones et al. 1120041 ). the "optimistic" scenario may vastly overestimate the true space density of 



interstellar comets. 



1.3. A New Approach 



The basic idea of this paper is to use white dwarfs with helium dominated atmospheres (DBs 
to limit the space density of interstellar comets. The large majority of whit e dwarfs have atmo- 



spheres which ar e hydrogen dominated and these stars are classified as DAs (jLiebert et al. 



1986 



Sion et al. 1 11988). The standard scenario is that gravitational settling of heavy elements is very 
effective, and therefore, with sufficient time, essentially all of the hy drogen is located in the s tar's 
outer envelope and dominates the composition of the photosphere (jHansen Sz Liebert 1 12003). A 
minority of white dwarfs have so little hydrogen that they have helium dominated atmospheres 
even though all their hydrogen resides in the star's outer envelope. If lines from neutral helium are 
detected, they are classified as DBs. If too cool for neutral helium to be detectable yet they display 
heavy elements in their spectra, they are classified as DZs. 

Two recent surveys show that white dwarfs accumulate hydrogen from an external source 
beca use there is a a correlation between t he average hydrogen mass and a DB's or DZ's cooling 



age (iDufour et al. 1120071 : IVoss et al. 1120071 ). The typical rate accretion rate is near 6 x 10 5 g s 1 



(jjura k, Xu Il2010l ). but there are large variations. In this pa per, we focus on t he ~50% of DBs 
with only upper bounds to the mass of accumulated hydrogen (IVoss et al. 1 120071 ). 



The mass of the convective zone where the accreted material resides is a strong function of 
the star's effective temperature. Because DBs with T* < 2 0,000 K have c onvective zones that are 
10,000 times more massive than DBs with T* > 22,000 K (jKoester Il2009l ). these two temperature 
classes may display external pollution quite differently. 

There is relatively little hydrogen in DBs warmer than 22,000 K. A conservative up per bound 

from optical and ultray iolet observations of most "warm" DBs is that n(H) /n(He) < 10 -4 ( Beauchamp et al. 



1999 



Voss et al. 



5 x 10 21 g (Koester 



2007). Since the mass of the outer convective zones of these warm DBs is at most 



20091 ). these stars have accumulated less than ~10 g of hydrogen. As de- 



scribed in Section 3, two particularly well studied DBs, GD 190 and BPM 17088, have accumulated 
less than 4 x 10 15 g of hydrogen, roughly the amount of this element carried in a single average-size 
comet. At the "typical" accretion rate of 6 x 10 5 g s _1 experienced by the cool DBs, these two 
stars would accumulate the upper bound to their hydrogen in less than 10 3 yr, much less time than 
their cooling ages of 50 Myr. Either the warm DBs do not accrete at the same rate as many cooler 
DBs, or they lose almost all their acquired hydrogen in a wind. 



- 4 - 



In Section 2 we elaborate in detail upon the "naive" idea that warm DBs have not accumulated 
much hydrogen simply because they have not accreted from ice-rich asteroids, comets or gas from 
the interstellar medium - the current best candidates as the sources for hydrogen in cooler DBs. 
While this argument hardly seems remarkable, there are interesting consequences. Using the toy 
model for a DB's hydrogen budget developed in Section 2, in Section 3 we constrain the mass carried 
within interstellar comets. In Section 4, we outline observational possibilities while in Section 5 we 
present our conclusions. 



2. TOY MODEL FOR A WARM DB'S HYDROGEN BUDGET 

2.1. Model Description 

We compute the amount of hydrogen in "warm" DBs for a regime when essentially all of its 
hydrogen of mass, M(H) has risen to the surface layers and is well mixed in the outer convective zone 
of mass M cvz . There may be a very thin non-convective photosphere lying above the convective 
zone, but we assume that the pho tosphere and the convective zone have the same composition 



(jDupuis et al. Hl992l : ISu h Lill2010l ). If hydrogen is accreted with rate M acc (H) and the total mass 



loss rate is M,,,,w, then: 



M(H) = M acc {H) - M wind (1) 



M cvz 



Equation (1) is mathematically equivalent to the expression governing the abundance of heavy 
atoms in the outer layer of a white dwarf where the atoms are "lost" because of downward settling 
(Koester 2009) instead of our case where hydrogen atoms are lost in an outward flow in a wind. 

For our purposes here, there is a crucial difference between accretion of gas from the interstellar 
medium and accretion from large solid objects such as asteroids. There is no reason to think that 
a quasi-spherical stellar win d strongly inhibits the accretion from a flat disk produced by tidal 



disruption of a parent body (I Jura M2003I ) . In contrast, such a wind can physically block interstellar 
gas accretion. In particular, while winds of 10 6 g s _1 can prevent interstellar accretion at the 
typical rate, once accreted hydrogen is diluted into a helium-dominated envelope, much higher 
wind outfows are required to inhibit the build-up of observable amounts of hydrogen. 

It is useful to introduce the characteristic time for the wind to deplete the mass in the convective 
zone, tdepi- 

tdepl = ■ (2) 

If M cvz , M w i n d and M acc are independent of time, the solution to Equation (1) is: 

M(H)[t] = M(H)[0] e-'/W + M acc (H)t depl (l - e" 1 ^') (3) 

Although Equation (3) provides a guide to our results, it is inexact because M cvz is not constant 
during the star's cooling. In the numerical solutions to Equation (1) presented below, we fit 
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the detailed calculati o ns fo r M cvz in a 0.6 M DB without hydrogen by iKoester I (|2009l ) and 



Gautschv fc Althaus I (|2002j ) in the temperature interval 22,000 K < % < 28,000 K with the 
expression: 



1.5 x 10 21 



25,000 



10 



(g) 



(4) 



Also, we relate the star's effective temperature to its cooling age by fitting the calculations of 
Bergeron^ for a 0.6 Mq DB white dwarf with: 



t « 49.7 



22, 000 V 



T, 



(Myr) 



(5) 



Because M cvz increases as the star ages, then, as can be seen from Equation (1), other factors being 
equal, the loss rate of hydrogen is diminished. Consequently, instead of reaching an asymptotic 
value as predicted by Equation (3), as illustrated in Figures 1-3, M(H) continues to increase with 
time. 



2.2. Boundary Condition 



We set t = for T* = 28,000 K, since at this moment in a DB's cooling, M cvz becomes large 
enough that Equation (1) can be sensibly employe d. The value of M(H)[0] is sensitive to the white 
dwarf's previous history (jUnglaub &; Bues 1 12000| ) . However, since the time for a 0.6 Mq DB to 
cool from T* = 28,000 K to T* = 22,000 K is about 35 Myr, then in most cases of interest, t > td ep i 
since we focus on wind loss rates greater than 10 6 g s _1 . If so, then the exact value of M(H)[0] 
is unimportant. In all cases, we adop t M(H)\0] = 1 x 10 15 g, a conservative upper bound for the 



K (ISullivan et al. 112008 ). 



ma ss of hydrogen in EC 20058-5234 (jPetitclerc et al. 112003 ). a well studied DB with % = 28,000 



2.3. White Dwarf Winds? 



By far, the largest uncertainty in Equation (1) is the stellar wind rate from white dwarfs. 



Unglaub I (|2008l ) has computed that radiatively-driven winds from white dwarfs with T* < 30,000 
K do not occur. Quite possibly, therefore, we should adopt M w i nc [ = 0. In contrast, in order to 



explai n the presence of small amounts of dredged-up carbon in warm DBs, iFontaine Brassard 
(2005) have proposed that M w i n d is as high as ~ 10 13 g s" 1 for T* > 20,000 K. However, the physical 



explanation for this putative wind is not provi ded. One possibility might be that the white dwarfs 



have hot coronae which could drive a wind (jLamers &; Cassinelli Ill999l ). In fact, white dwarfs 



do have convective envelopes and there could be some transfer of the convective energy flux into 
heating the outer layers of the gas. 



1 see http://www.astro.umontreal.ca/~bergeron/CoolingModcls/Tables/Table_DB 
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We use two lines of evi dence to bound the w ind rates from warm white dwarfs. First, these 
stars are not X-ray sources (jO'Dwver et al. II2003I ). Assu me a wind outflow speed n ear the white 
dwarf's characteristic escape velocity, V esc , of 4000 km s" 1 ( Lamers Cassinelli 1 ll999l ) and consider 
the kinetic energy flow carried by the wind, 1/2 M w i n d V^ sc . Much of this energy is deposited into 
the surr ounding interstellar me dium to create a hot bubble. Using the result that Lx < 4x 10 



erg s" 1 (IMusielak et al. 1 120031 ) for GD 358, a DB with T* = 27,000 (IProvencal et al. 1120091 ). it 



seems that M w i nc i <~5x 10 g s . (This limit could be relaxed if the X- r ays ar e deposited in 
a region greater than the 5'.'1 circle on the sky employed by IMusielak et al. I (|2003l ) in their data 
analysis. However, as can be seen from their Figure 2, there is no extended X-ray emission.) 

Second, some white dwarfs accrete and therefore are not losing a large amount of mass in a 
wind. Consider SDSS 1228+1040 with T* = 22,000 K, the hottest white dwarf known to be heavily 
polluted and possess a dust disk (Gaensicke et al. 200 6, Brinkworth et al. 2009, Melis et al. 2010). 
Because this DA does not have a convective envelope (jKoester Il2009l ). Equation (1) does not apply 
Instead, to accumulate heavy elements, it is necessary that M w i n d < M acc . From the formalism 



and p arameters given by iKoester I (|2009l ). the magnesium abundance reported by lGaensicke et al. 



(200 7]) and the pres umption that ~10% of the mass of the material is magnesium (see, for example, 
Klein et al. I ( 2010 )). then the heavy atom accretion rate is about 2 x 10 9 g s _1 . This value is 



within the range of 3 x 10 8 g s 1 to ~ 10 1 g s 1 infe rred for most white dwarfs with dust disks 



(jKilic et al. 1120061 ; Ivon Hippel et al. 



2007 



Farihi et al. Il2009l ). It seems likely that the wind rate 



from this star is less than 10 10 g s _1 . (We assume that this DA without a convective envelope has 
the same mass loss rate as a DB with a low mass convective envelope. Therefore, this argument is 
not completely certain.) 

The source of the hydrogen in DBs with T* < 20,000 K Is uncertain. If some o f these stars 



accret e interstellar gas with the typical time- averaged hydrogen accretion rate given bv lVoss et al. 
(120071 ) of ~10 6 g s _1 , then M W i n d likely is less than this value. Because warm DBs have much less 
massive convective envelopes than cool DBs, there is no reason to imagine that M W i n d increases 
dramatically with effective temperature. We conclude that the wind rates from warm DBs likely 
are less than 10 10 g s _1 and plausibly less than 10 6 g s _1 . 

Besides a wind, white dwarfs also may lose mass in explosive events as occurs in novae. Because 
there is no evidence that isolated white dwarfs undergo such events nor is there any theoretical 
reason to expect them, we ignore this possibility. 



2.4. Interstellar Gas Accretion? 



Although considered for many years, the physics of accreting interstellar hydrogen gas by 
white dwarfs is not well understood. It seems, for example, that accretion at the Bondi-Hoyle 



rate is excluded (jFarihi et al. II2010I ). DBs younger than 50 Myr have performed less than one full 



oscillation through the Galactic Plane. Therefore, these stars may have only encountered hot, low 
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density material and consequently have accreted interstellar gas at a negligible rate during their 
entire cooling age. Because accretion of interstellar gas-phase hydrogen can only strengthen the 
upper limits to asteroidal and cometary accretion discussed below, we neglect this possibility. 



2.5. Accretion of asteroids with ice? 



In our calculations for asteroidal accretion, we assume that M acc (H) = 6 x 10 5 g s" 1 , a 
characteristi c hydrogen accreti on rate for DBs with T* < 20,000 K although there is a large variation 
in this rate (|Voss et al. 1 120071 ). We show in Figure 1 more detailed solutions to Equation (1). In 
these models, we see that if M win( i < 10 10 g s _1 , then DBs typically accumulate more hydrogen 
than observed. This means that DBs warmer than 22,000 K either have "strong" winds or they do 
not accrete from ice-rich asteroids. 



2.6. Oort Cloud Analog Cometary Accretion? 



Alcock et al. I (119861 ) proposed that direct impact by comets can produce observable signatures 
in a white dwarf's photosphere. Here, we estimate the hydrogen p ollution of a D B white dwarf by 
comets in an analog to the solar system's Oort cloud. According to lFrancis I (120051 ). 0.8 dynamically 
new comets AU _1 yr _1 arrive in the inner solar system. With a flat distribution of perihelia, the 
expected rate of collisions onto a star of radius 0.013 R0 is 4.8 x 10 -5 yr" 1 . With an average 
comet mass somewhere between 5.6 x 10 16 g and 1.2 x 10 18 g (Francis 2005) and the assumption 
that 5% of the comet's mass is hydrogen (jGreenberg Ill998l ). then the hydrogen accretion rate is 
expected to lie between 4 x 10 3 g s _1 and 9 x 10 4 g s . To be conservative, we adopt the lower 
bound and show results from Equation (1) in Figure 2. We see that if M w i n d < 10 8 g s _1 , then 
as long as T* > 22,000 K, the star accumulates more hydrogen than allo wed for the typical upper 
bound for DBs obtained from optical spectra (jBeauchamp et al. Ill999l ), Therefore, if the winds 
from these white dwarfs are plausibly weak, the mass of their analogs to Oort cluds are smaller 
than the solar system's. Either these stars never had many pri moridal comets or they lo st their 
comets when the star evolved on the Asymptotic Giant Branch (IParrtiott &: Alcock Ill998l ). 



We assume that the comet's hydrogen is visible in the star's outer convective zone. In fact, 
during its high-spee d impact, the iceball might penetrate beneath the convective zone. Following 
Alcock et al. I (| 19861 ) , we assume that the material vaporized from the comet has such a high entropy 
that it quickly rises to the surface of the star and therefore could be detected. 
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3. LIMITS ON INTERSTELLAR COMETS 

Having argued above that the well known sources of hydrogen accretion are not important in 
many warm DBs, we now proceed to estimate the space density of interstellar comets. 



3.1. Computed Accretion Rate 

We assume that the comets have hyperbolic orbits with speed relative to the white dwarf of 
V w d- If the comet impacts a white dwarf of radius and mass M*, then it is 100% accreted; 
we assume there is no accretion from a near-miss. We treat the ambient interstellar iceballs as a 
continuous fluid with a hydrogen mass density, pb a u(H). Because the escape speed from the surface 
of the white dwar f is much greate r than V w d, then the total mass of accreted hydrogen, M acc (H) is 
(see, for example. ISpitzer I ()1978l )): 



M acc {H) 



Pball(H)[t] 

V wd (t) 



dt 



(6) 



Below, we describe how pb a u and V w d may vary with time. 



Even before impact, incoming comets may split as they a re heated and /o r tidally disrupted. 
However, the resulting fragmentation speeds of ~0.1 km s _1 ( Sekanina 1 1 19821 } are vastly smaller 
than the impact speeds that exceed 1000 km s . We therefore neglect the small deviations from 
the initial hyperbolic orbits that arise. In any case, the loss of some fragments that deviate away 
from a path that leads them into a direct collision is balanced by the gain of other fragments from 
parent bodies on a near-miss trajectory that deviate into a collision-course. 

For simplicity, we assume V w d is dominated by the white dwarf's space motion in the Local 
Standard of Rest and we ignore the motions of the comets. Assuming that only the motion vertical 
to the plane varies significantly with time, we write: 



v wd = {U 2 + V 2 + W(t) 



2\V2 



(7) 



with U, V and W having their conventional meanings. We adopt the simplify ing assumption that 
relat ive to the Galactic Plane a star undergoes vertical simple harmonic motion (jKuiiken Gilmore 
1989h 

W(t) = W cos k 1 ' 2 1 (8) 



With k 1/2 = 2.3 x 10 -15 s" 1 (jSiebert et al. 112003 ). the period of a star's oscillation through the 
Galactic Plane is 86 Myr. 

We assume that interstellar iceballs have the same spatial distribution as F-type stars and 
therefore relative to the Galactic Plane: 



Pbaii{H)[z] = Pball (H)[0]e 



-\z\/h 



(9) 
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We adopt h = 150 pc (|Gilmore Zeilik Il2000h 



Our analysis only applies to iceballs in a limited size range. Small iceballs are significantly 
sublimated on their approach to the white dwarf, an d the result ing gas atoms may not accrete 
onto the star. Reformulating slightly Equation (14) in Jura (j2005l ). the thickness, AR, of a comet 
which is thermally sublimated as it approaches a star is: 

A R ~ L * m H2Q^subl ( m 

~ 167rAS[2 G R* M*] V 2 { ' 

where mn 2 o is the mass of a water molecule, A<j) su yi is the angular variation in its orbit traced by 
the incoming comet as seen by the host star where ice sublimation dominates the cooling of the 
whit e dwarf, AE is t he energy per molecule associated with sublimation (taken equal to 2 x 10~ 12 
erg. ISekanina I ([20020 ) and we assume an orbit that is nearly parabolic. We take Acp su bi = n/2 since 
ice sublimates significantly at distances much greater than i?*. For illustrative purposes, we take 
L* = 0.05 L Q , i?* = 0.013 R and M* = 0.6 M Q as appropriate for the white dwarfs of interest. 
We find that AR ~ 240 cm. Therefore, iceballs with radii greater than ~10 m deposit most of their 
mass onto the white dwarf. Our analysis only constrains iceballs with hydrogen masses less than 
the upper limits of accreted hydrogen in our target white dwarfs. This implies that we consider 
iceballs with less than 4 x 10 15 g of hydrogen and therefore a total mass less than 8 x 10 16 g and 
a radius less than 2 km. 

Besides thermal sublimation, is there some other process which leads to a comet's disintegration 
and would the resulting gas and debris somehow be inhibited from accreting onto the star? In 
considering this possibility, it should be recognized that the approach of a comet to a white dwarf 
is less disruptive than the approach of a comet to the Sun for two reasons. First, the white dwarf's 
luminosity is lower than the Sun's and thus a comet nearing a white dwarf is heated less and the 
resulting thermal stresses are smaller. Second, the typical white dwarf mass of 0.6 M is less than 
the Sun's mass and therefore the gravitational tidal forces from a white dwarf are smaller. As a 
result, while comets may split as they approach their target white dwarf, analogous to the approach 
of comets to the Sun, we do not expect them to fully disintegrate - if they ever do - until they 
are closer than ~ 1 Rq. At this separation, the comet's incoming speed is nearly 500 km s _1 and 
it survives less than 1000 s before impact. Perhaps as the comet nears within 1 Rq of the white 
dwarf, there may be extensive explosions. Generaliz ing from the observ ations of comet 17P/Holmes 



which displayed remarkable transient brightening (jReach et al. 1120100 . the typical outflow speed 
of dust fragments and gas debris in such events might be 0.1 km s _1 . At this speed, the ejected 
material travels ~100 km in the comet's frame of reference before impacting the white dwarf. 
Since this transverse distance is much smaller than the radius of the white dwarf of nearly 10 4 
km, we expect that cometary material that initially is on a trajectory aimed at the star does in 
fact impact, regardless of any explosions that occur during the last 1000 s. This argument applies 
even to small dust grains. The luminosity of a white dwarf is so low tha t even tiny solid particl es 



cannot be substantially deflected by radiation pressure (see Section 5.6 of lPlavchan et al. I (|2005l )) 



In conclusion, even if a comet largely disintegrates during its final approach to a white dwarf, the 
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bulk of the material is accreted. 



3.2. Application to GD 190 and BPM 17088 



To date, the most sensitive upper limits to hydrogen in DBs are by iPetitclerc et al. I (120051 ) 
who report ultraviolet spectroscopy for three stars. One of the stars, EC 20058-5234 with T* 
28,000 K, is relatively y oung and just beginning to have a well developed convective zone 



(jGautschv Althaus 1120021 ), Our toy model described in Section 2 may not apply. Therefore, 
we only quantitatively discuss results for GD 190 and BPM 17088. The parameters that we adopt 
for the two targets and relevant references are given in Table 1. If their effective temperatures are 
less than 22,000 K, then our derived upper limits to the space density of interstellar comets are 
weakened. 

From Equation (8), the maximum distance that a star travels from the Galactic Plane is 
Wo/k 1 ? 2 . If W equals the value of W given in Table 1, then GD 190 and BPM 17088 remain 
within 150 pc, the estimated scale height around the Galactic Plane of the interstellar iceballs. 
Using the parameters in Table 1, and making the simplifying assumption that the interstellar 
density and speed of the white dwarf are constant since |W(i)| << V w d and \z\ < h, then from 
Equation (6), the upper limits to Pbaii(H) experienced by GD 190 and BPM 17088 are 1.5 x 10 -29 g 
cm -3 and 1.0 x 10 -29 g cm -3 , respectively. We assume that interstellar comets are mainly ice and 
silicates and that p(0) 10 p(H) (Greenberg 1998). Consequently the upper bounds to PballiP) 
in the Galactic midplane are 1.5 x 10~ 28 g cm -3 and 1.0 x 10 -28 g cm -3 for GD 190 and BPM 
17088, respectively. Therefore, we find that the density of oxgyen contained within large iceballs 
is less than about 1% of the total interstellar oxygen density, and that "optimistic" models for the 
space density of interstellar comets (Section 1.2) with radii between 10 m and 2 km are excluded. 

Above, we assumed that any wind from the white dwarf is completely negligible. We now 
consider the implications of including a wind and show in Figure 3 results for the amount of 
hydrogen accumulated by a DB as a function of time for a case where the average accretion rate 
is 10 g s _1 , a factor of four above the values given in Table 1. We see that if M W i n d is less than 
10 6 s _1 , more hydrogen is accumulated than allowed by the observed upper bound of 4 x 10 15 g. 
Therefore, if warm white dwarf winds are plausibly weak, interstellar iceballs are rare. 



4. OBSERVATIONAL POSSIBILITIES 



We have suggested that DBs warmer than 22,000 K are not polluted by ice-rich asteroids, 
comets or interstellar gas. This s cenario can be tes t ed by examining atmospheric compositions. 
Using ground-based observations, IZuckerman et al. I (|201Q ) found that over 30% of DBs between 



13,500 K and 19,500 K display heavy atom s. To date, the most extensive similar study of DBs 



with T* > 22,000 K is bv IVoss et al. I (|2007h . In this sample of 11 stars, only WD 2354+1659 was 
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identified as having heavy atoms in its photosphere. One possibility is that this star has accreted ice- 
poor asteroids as c ould be inferred with a detailed abundance analysis similar to that performed by 



Klein et al. I (|2010l ) for GD 40. Also, however, the temperature of this star is uncertain. Depe nding 



upon the gravity and hydrogen conten t, T* could be as high as 24,800 K (|Voss et al. 1 120071 ) or as 



low as 19,000 K ([Koester et al. 1120051 ). While the frequency of externally-polluted "warm" DBs 



appears to be low comp ared to the c ooler DBs, further measurements should be performed since 
th e spectra obtained bvlVqss et al. I (j2007l ) are not as sensitive to heavy atoms as those obtained 



bv lZuckerman et al. I (|2010l ) 



Ultravio let spectroscopy is a se nsitive tool to measure heavy elements in a white dwarf's 



atmosphere (jDesharnais et al. 1 12008). Data from the Far Ultraviolet Spectroscopic Explorer and 



the Hubble Space Teles cope have been obtain ed for 5 DBs wa rmer than 22,000 K: G D 190, BPM 
17088. EG 20058-5234 (Fetitclerc et abl hoosh. PG 0112+104 dProvencal et al. Ihoooh and GD 358 



([Provencal et al. II1996I ). Although these stars possess atmospheric carbon, it is almost certain that 



this heavy element is intrinsic to t he star and not accreted (Fontaine & Brassard 2005, Desharnais 
et al. 2008). With t he mo dels of iKoester I (120091 ) . the upper limits to the iron abundance found 



by iPetitclerc et al. I (|2005l ) and th e assumption that 1/3 of the accreted mass is iron as seems 



appropriate for extrasolar asteroids ([Klein et al. IZuckerman et al. I l2007l . l20ld : bufour et al" 



2010 



Vennes et al. Il2010l ). then in GD 190, for example, the upper limit to the h e avy at om accretion 



rate is ~1 x 10 6 g s , a value lower than any reported bv lZuckerman et al. I (|2010l ) for polluted 
cooler DBs. Therefore ultraviolet results are more sensitive than optical studies for identifying 
pollution. If each ultraviolet targeted DB has at least a 30% probability of experiencing pollution, 
as found by IZuckerman et al. I (|2010l ) for cool DBs, then the total probability that all 5 atars are 
uncontaminated is less than (0.7) 5 or 16%. Future observations of warm DBs may further test the 
prediction that these stars do not accrete ice-rich asteroids. 



The warmest DA known to have a dust disk i s SDSS 1228+1040 with T* = 22,600 K ([Gaensicke et al 



2006 



Brinkworth et al. 



2009 



Melis et al. Il2010l ) . Discovery of additional warm DAs with dust and 



then measuring their pollution may help extend the argument described in Section 2 to constrain 
the mass loss rates from such stars. 

In our toy model described in Section 2, we have treated hydrogen accretion as a minor 
perturbation in a DB atmosphere. In fact, it may take only ~ 10 18 g of hydrogen to prevent the 
conversion of a DA into a DB for stars wa rmer than 22,000 K by dredge-up of interior helium 



( MacDonald fc Vennes 



1991 



Su fc Li Il2010l ). In contrast, for cooler stars with much more massive 



helium convective zones, considerably more hydrogen may be required to prevent the conversion 



from DAs to DBs ([MacDonald & Vennes Ill99ll : iTremblay et al. 1120081 ). The fraction of all white 



dwarfs that are DBs with T* between 25,000 K and 40,000 K is a factor of 2.5 les s than the same 
fraction of DBs with % between 16,000 K and 22,000 K (jEisenstein et al. II2OO6I ). We speculate 
that the systems with ice-rich asteroids accrete enough hydrogen to prevent a white dwarf from 
becoming a DB until the temperature falls appreciably below 22,000 K. 
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Our supp osition that th e Por t cloud is unusually massive is consistent with the previous 
suggestion by IShull k Stern I (|1995l ). and the picture that the solar system is som ewhat special. 
The orbital eccentricities of our planets are unusually low (IMarcv k Butler 2000J) and extrasolar 



analogs to the Late Heavy Bombardment appear to be rare (jBooth et al. 1120091 ) 



5. CONCLUSIONS 

We identify a subclass of warm DBs where accretion of asteroids, analogs to Oort clouds comets 
and interstellar gas all appear to be negligible. We make the plausible argument that the wind mass 
loss rate these stars is less than 10 6 g s . If these stars have their nominal effective temperatures, it 
then follows that because they have accumulated only a small mass of hydrogen, interstellar comets 
with radii between 10 m and 2 km a contain less than 1% of all interstellar oxygen. "Optimistic" 
scenarios for the formation of the solar system's Oort cloud may not apply to most other stars. 

This work has been partly supported by the National Science Foundation. I thank B. Klein, 
M. Kilic, D. Koester, S. Xu and B. Zuckerman for helpful comments. 
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Table 1 - Target White Dwarfs 





GD 190 


T 1 T 1 A if 1 prrinO 

BPM 17088 


Rererence(s) 


rri /t/'\ 

T* (K) 


22,000 


oo r* r\r\ 

22,600 


IX 


log g (cgs) 


7.97 


8.07 


ab 


"71 T / ~\ IT \ 

(M ) 


0.59 


0.64 


ab 


1 /TT\ / / T T \ 

log n(H)/n(He) 


< -5.5 


< -5.5 




log ikZcvz (g) 


Zl.io 


21.01 


c 


log M occ (iT) (g) 


<15.63 


< 15.51 




log i coo « (Myr) 


7.70 


7.64 


d 


logi?* (cm) 


8.96 


8.93 




M acc (H) (g s^ 1 ) 


<2.7 


<2.4 




(km s" 1 ) 


-19 


14 


e 


V (km s" 1 ) 


-12 


-10 


e 


W (km s" 1 ) 


-9 


8 


e 


V wd (km s" 1 ) 


24 


19 





" IPetitclerc et al. I (|2005h . Woss et al. I (120071 ). c Equation (4) and ignoring t he mass var i ation of this 
quantity, ^Equation (5) and ignoring the mass variation in this quantity, iSion et al. I (1988) 
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T. (1000 K) 

Fig. 1. — Predicted accumulated hydrogen mass from asteroidal accretion as a function of effective 
temperature for different wind mass loss rates. We use M(H)[0] = 1.0 x 10 15 g, M acc (H) = 6 x 
10 5 g s _1 and values of M w i n d of 10 10 , 10 9 and 10 8 g s -1 shown by the red, blue and green lines, 
respectively. The dotted line represents the upper limit to hydrogen in a DB with the assumptions 
that n(H)/n(He) < 10 -4 and M cvz given by Equation (4). We see that in the likely case that M W i n d 
<10 10 g s _1 , the assumption that the star accretes ice-rich asteroids over-predicts the amount of 
accumulated hydrogen 
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Fig. 2. — Similar to Figure 1 except for accretion from an analog to the Oort cloud with M acc (H) 
= 4 x 10 3 g s _1 and values of M W i n d of 10 8 , 10 7 and 10 6 g s -1 shown by the red, blue and green 
curves, respectively. We see that in the plausible case that M w i n d < 10 8 g s _1 , the assumption that 
the gas accretes comets from a conservative extrapolation of our own Oort cloud over-predicts the 
amount of accumulated hydrogen. 
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T. (1000 K) 

Fig. 3. — Similar to Figure 1 except for an accretion rate of M acc = 10 g s^ 1 and values of M w ind 
of 10 7 , 10 6 and 10 5 g s _1 shown by the red, blue and green curves, respectively. The downward 
pointing arrow represents the upper bound to the mass of accumulated hydrogen in GD 190 at the 
temperature range of the two DBs in Table 1. We see that in the plausible case that M w ind < 10 6 
g s _1 , more hydrogen is predicted to be accumulated than observed. 



